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ABSTRACT 



We investigate molecular evolution from a molecular cloud core to a first hy- 
drostatic core in three spatial dimensions. We perform a radiation hydrodynamic 
simulation in order to trace fluid parcels, in which molecular evolution is investi- 
gated, using a gas-phase and grain-surface chemical reaction network. We derive 
spatial distributions of molecular abundances and column densities in the core 
harboring the first core. We find that the total of gas and ice abundances of many 
species in a cold era (10 K) remain unaltered until the temperature reaches ~500 
K. The gas abundances in the warm envelope and the outer layer of the first core 
(T < 500 K) are mainly determined via the sublimation of ice-mantle species. 
Above 500 K, the abundant molecules, such as H 2 CO, start to be destroyed, and 
simple molecules, such as CO, H 2 and N 2 are reformed. On the other hand, 
some molecules are effectively formed at high temperature; carbon-chains, such 
as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also 
find that large organic molecules, such as CH3OH and HCOOCH3, are associated 
with the first core (r < 10 AU). Although the abundances of these molecules in 
the first core stage are comparable or less than in the protostellar stage (hot 
corino), reflecting the lower luminosity of the central object, their column densi- 
ties in our model are comparable to the observed values toward the prototypical 
hot corino, IRAS 16293-2422. We propose that these large organic molecules can 
be good tracers of the first cores. 

Subject headings: Astrochemistry — ISM: molecules — ISM: clouds — Stars: 
formation 
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INTRODUCTION 



It is well established that a star is formed by gravitational collapse of a molecular 
cloud core. The collapsing core is initially optically thin to the thermal emission of dust 
grains, and undergo isothermal run-away collapse as long as the cooling rate overwhelms 
the compressional heating. The isothermal condition breaks down when the central 
density reaches ~10 10 cm -3 , and the temperature starts rising. Increasing gas pressure 
decelerates the contraction, and th e core come to the hydrostatic equilib rium, which is 
called a first hydrostatic core (e.g., lLarson Ill969l ; iMasunaga et al. Ill998l ). Although the 
first core is a transient object, it plays essential roles in determining the phys ical condition 
of a protostar; it fragments and forms binar ies (IMatsumoto &: Hanawa 1120031 ). and drives 
bipolar molecular outflows (ITomisaka 1 120021 ) . When the temperature reaches ~2000 K and 
hydrogen molecules start to dissociate, the central region of the first core collapses again to 
form the protostar. 

Star formation processes include both the structure evolution as mentioned above, and 
molecular evolution. Rates of chemical reactions depend on various physical quantities. 
Among them, the temperature is crucial; e.g., rates of ice sublimation and reactions with 
potential barriers are proportional to exp(— AE/kT), where AE represents desorption 
energy or a height of an energy barrier. In addition, the molecular evolution in star forming 
regions is a non-equilibrium process. Therefore, accurate temperature determination via 
radiation hydrodynamic (RHD) simulations is important to investigate molecular evolution 
in star forming cores. 

There are a few previous works which combined a chemical reaction network model 
with a RHD model of star-forming cores. Aikawa et al. (2008, here after AW08) 
investigated chemistry from a mol ecular cloud core to a protos tellar core adopting a 
spherically symmetric RHD model ( IMasunaga Sz Inutsuka 1120001 ). and mainly discussed the 
molecular evolution that occurs at r > 10 AU, where T < 300 K. The spherical symmetry, 
however, should br ake down inside the centr ifugal radius of ~100 AU, where circumstellar 
disks will appear. Ivan Weeren et al. I (120091 ) adopted an axis symmetric RHD model 



( Yorke fc Bodenheimer 19991 ). and investigated molecular evolution mainly at T < 100 



K from a collapsing molecular cloud cores to a forming circumstellar disk. The spatial 
resolution of their physical model is ~3 AU in the central region. So the central first core 
and/or protostar is not resolved. 

In this paper, we combine a detailed chemical reaction network model of gas-phase 
and grain-surface chemistry with a three dimensional RHD simulation of a star-forming 
core for the first time, and investigate molecular evolution at T = 10-2000 K. So far, a few 
groups have successfully performed three-dimensional R(M)HD simulations of a low mass 
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star- forming core up to the first core stage (jWhitehouse fc Bate 1120061 ; ICommercon et al. 



20101 ; iTomida et al. H2010al Jbl). On the other hand, evolution beyond first c o res in three 



dimension is m ore c hallenging and investig ated only very recently by iBate I (120101 . l201ll ) 



Tomida I (120 121 ). and ITomida et al. I ( 120121 ). So we focus on chemistry in the first core stage 



in the present work. Our motivations are twofold. 

Firstly, first cores are good targets for Atacama Large Millime ter /submillimeter Array 
(ALMA). Although over 40 years have passed since lLarson I fjl969h predicted the presence 
of first cores, they have not yet been detected. Observation of first cores is challenging; 
they are buried in dense envelopes, and have very compact structure (~10 AU) and short 
l ifetime (a few 100 yr). Recentl y, some candidates have been report ed: IRS2E in L1448 



(IChen et al. Il2010[ ). Per-Bolo 58 ( Enoch et al. Il2010h and L1541-mm (IPineda et al. 1 1201 if ) 



Observational properties of these objects, low bolometric luminosity L < 0.1 Lq and cold 
spect ral energy distribution, are consistent with the t heoretical predictions of first cores 



[e.g. 



Masunaga et al. 



19981 : ISaigo fc Tomisaka 11201 lb . On the other hand, IRS2E has 
a high velocity outflow (~25 km/s), and Per-B olo 58 is considerably luminous at 24 /xm 



and associated with a well-collimated outflow ( iDunham et al. 1 1201 if ), which contradict 



the theoretical prediction that the ou tflow from first cores is slow (<5 km/s) and not 



well collimated ( IMachida et al. 1 120081 ). Confirmation and further investigations of the 



candidates require observation of molecular emission lines. Therefore, chemical models of 
first cores are highly desired. 

Secondly, recent hydrodynamic simulations have shown that the outer region of first 
cores might dire ctly evolve to circumstellar disks, while the central region collaps es to 
form protostars JSaigo et al. I I2OO8I : IMachida et al~ll201o[ l201l[ Bate I boiol . boilh . In 
that case, the chemical composition of the first cores would be the initial composition of 
the circumstellar disks. The Molecular evolution during the formation and evolution of 
the circumstellar disks is still an open question. Although most chemical models of the 
circumstellar disks have used abundances of the molecular cloud cores as initial abundances, 



its validity is questioned (Ivan Weeren et al. 1120091 ; IVisser et al. 1 1201 if ). In other words, 
there remains a missing link between chemistry in the molecular cloud cores and the 
circumstellar disks. Our detailed three-dimensional chemical-hydrodynamical model from 
the molecular cloud cores to the first cores is the first step to fill the gap. 

The rest of the paper is organized as follows. In Section 2, we briefly describe our RHD 
simulations and chemical reaction network model. In Section 3, we discuss trajectories of 
fluid parcels in our RHD simulations. We show molecular evolutions in the assorted fluid 
parcels, and the spatial distributions of molecular abundances in the first core stage. In 
Section 4, we discuss molecular column densities and uncertainties in our reaction network 



-5 - 



model. We compare our model with models of hot corinos in terms of physical and chemical 
properties. We also discuss the effect of accretion shocks on chemistry, and chemistry after 
the first core stage. We summarize our results in Section 5. 



2. MODEL 

2.1. Physical Model: A 3 Dimensional Radiation Hydrodynamic Simulation 

Physical structure of a collapsing molecular cloud core is calculated using a three- 
dimensional neste d-grid self-gray i tationa l radiation magnetohydrodynamics code, which 
was developed by Tomida et al. (j2010a ). Here we briefly describe physical parameters and 



structure of a core. 

The self-gravitational magnetohydrodynamics equations are coupled with ra- 
diation transfer, which is treat ed with gray flux-limited diffusion approximation 



( iLevermore fc Pomraning Ill98ll ). The gas temperature is set to be equal to the dust 
temperature. In the present work, the magnetic field is set to zero for simplicity. Idealized 
equation of state is assumed with the adiabatic index of 5/3. The simulation consists of 
multiple nested-grids in order to achieve higher resolution in the central region. Each level 
of the nested grids has 64 3 cubic cells in (x, y, z). At the end of the simulation, 13 levels of 
nested-grids are generated and the smallest scale corresponds to ~0.05 AU. 

We calculate the evolution of a molecular cloud core with the mass of 1 M Q . The 
initial density distribution is given by a critical Bonnor-Ebert sphere, and enhanced by 
a factor of 1.6 so that the sphere becomes gravitationally unstable. Initially, the central 
number density of the core is n c ~ 8.3 x 10 5 cm -3 (p c = 3.2 x 10~ 18 g cm" 3 ). The radius 
of the cloud core is 6300 AU. The initial angular velocity is 2.0 x 10~ 14 s _1 , and the initial 
temperature is set to 10 K. The simulation is stopped when the maximum temperature 
reaches 2000 K, at which temperature hydrogen molecules collisionally dissociate, and the 
second collapse starts. 

In this paper, we define the moment when the central gas number density of the core 
reaches n c ~ 2.6 x 10 10 cm -3 (p c = 10~ 13 g cm -3 ) as tpc = 0; the first core is born. Figure 
[1] shows the temporal variations of the central gas number density and temperature of 
the core at tFc > 0. The temperature and the density increase by two and five orders of 
magnitude in ~3000 yr, respectively, which is the lifetime of the first core in our model. 
Figure |5] shows the core structures at tpc = 1500 yr (middle stage) and 2800 yr (late 
stage). We present molecular distributions in these two stages in Section 3. The core is 
oblate or disk-like because of the rotation. In the panels (c) and (f), the dashed lines depict 
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temperature contours of 25 K and 100 K, which are roughly the sublimation temperature of 
CO and large organic molecules, respectively. As the core evolves from 1500 yr to 2800 yr, 
the sublimation radii of CO and large organic molecules expand from ~40 AU to ~100 AU 
and from ~5 AU to ~10 AU, respectively. We discuss the physical structures of the core in 
more detail in Section 3.1. 

Chemical processes are not included in the hydrodynamic simulation. In this work 
we calculate molecular abundances as a post process; we trace 10 5 fluid parcels in the 
hydrodynamic simulation, and calculate a chemical reaction network model along 10 4 



selected trajectories. Our approach is similar to AW08 and Ivan Weeren et al. I (120091 ) . 
Initially, the fluid parcels are distributed spherically at r = 1000-2500 AU. The parcels 
follow the local flow at each hydrodynamic time step. Leapfrog integration is used for 
tracing the parcels, and the physical parameters at the position of the parcels are obtained 
by the trilinear interpolation. In the calculations of molecular evolution in the fluid parcels, 
the temperature, density, visual extinction, and cosmic-ray ionization rate are updated at 
each time step along the trajectories (see also Section 2.2.4). We note that the accretion 
shock heating layer at the first cor e surface is not spatially resolved in our model, since the 
thickness of the layer, ~10~ 3 AU (ISaigo &: Tomisaka Il201ll ). is smaller than the smallest 
scale of our calculation. Effect of the shock layer on chemistry is discussed in Section 4.4. 



2.2. Chemical Reaction Network 



The dynamical evolution covers a wide range of temperature and density: 10-2000 K 
and 10 5 -10^~ 

"* fl2006h 



condit ions, we com b ine tw o chemical reaction network models, iGarrod fc Herbst 



and lHarada et al. I (120 10[ ). as a base of our model. The former is a gas-grain reaction 
network model applied to a star-forming core at 10-200 K. The latter is a high temperature 
gas-phase reaction network model, and originally developed for the temperature of 100-800 
K. The major difference between the two network models is that reactions with a high 
potential energy barrier are included in the latter model. For the neutral-neutral reactions 
and neutral-ion reactions, we switch the rates from the former network to the latter network 
when the temperature exceeds 100 K. In the following, we describe modifications to gas 
phase reactions, treatment of gas-grain interactions, and grain-surface reactions in our 
model. In total, the reaction network consists of 461 gas-phase species (including dust 
grains with three different charge states) and 195 grain-surface species, an d 11696 reaction s. 



The numerical code for sol ving; rate equations is adapted from Nautilus ( IHersant et al 



20091 ; ISemenov et al. 1120101 ) 
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2.2.1. Modifications to Gas Phase Reactions 



Firstly, we modi fy branching ratios of dissociative recombination reactions. 
Geppert et al. I (120061 ) suggested that formation of saturated complex molecules by 
dissociative recombination of their protonated precursors is not efficient. For example, 
when protonated methanol (CHsOHj") rec ombines, the branching ratio for CH 3 OH + H is 
only 3%. Recently, lHamberg et al. I ( 120101 ) also showed that in dissociative recombination 



of dimethyl ether cation (CHsOCH^J"), the branchin g ratio for CH3OCH3 + H is only 



7%. Referri ng to these experiments, the network of iGarrod et al. I (120081 ). which is an 



extension of IGarrod &: Herbst I (120061 ) to more complex species, assumed that two- fragment 
channels (ex. X + + e~ — > Y + Z) represent 5% each, and that the remai nder is evenly split 



betwe en channels with three or more fragments, while in the network of IGarrod fc Herbst 



(120061 ) the branching ratios of dissociative reactions of protonated complex molecules are 
equally weighted for most species. We modify the branching ra tios in the recombinat ion of 
protonated precursors of saturated complex molecules following IGarrod et al. I (120081 ) . 



We add some three-body association reactio ns and collisional diss ociation reactions 



for all neutral species to our network referring to IWillacy et al. I ( 



998 ) (see App e ndix A ) . 



Harada et al. I (120 lOf ) 



Collisional dissociation of H2 is an exception; we use the rates in 

instead of I Willacy et al. I (119981 ). Since collisional dissociation is endothermic, the rate of 
the collisional dissociation reactions is proportional to exp(— 7/T). The energy barrier, 
7, must be larger than AH, which is the dif f erence between the formation enthalpy 
of products and reactants. In IWillacy et al. I (119981 ) the value of 7 is set to be lower 



than AH for some reactions. For those reactions, we calculate AH at 300 K referring 



to NIST-JANAF Thermochemical Tables ( |http: / /kinetics.nist.gov/janaf/ ) and Kinetic 
Database for Astrochemistry (KIDA; http://kida.obs.u-bordeauxl.fr; Wakelam et al. 
2012), and rep l ace 7 with AH. For the species which are not listed in the network of 
Willacy et al. I (119981 ) the collisional dissociation rates are set to be, 



k cd = 1.0 x 10- 10 exp(-A#/T). (1) 

If the formation enthalpy of reactants or products are not available in the literature, we 
assume AH is 5 x 10 4 K (~5 eV). We assume that the collisional dissociation reactions 
preserve functional groups. Otherwise, we assume that one H atom is extracted via the 
collisional dissociation. It should be noted that the set of collisional dissociation reactions 
in our model is likely incomplete because of the lack of experimental data; for many 
reactions listed in Table A2, there is no laboratory experiments in which products are 
confirmed. If the products are different from those listed in Table A2, the value of 7 should 
also be different. In Section 4.2, we discuss how the molecular abundances depend on the 
assumption on the collisional dissociations. 



-8 - 



2.2.2. Gas- Grain Interactions 



In molecular cloud cores, dust grains are negatively charged. But the p opulation 
of charged particles changes as the gas density increases (jUmebayashi 1 1 1 9 8 3h . Grains 
become the dominant charged particles at high densities (n > 10 10 cm -3 ), and eventually 
grain-surface recombination (ion-grain collision) overwhelms gas-phase recombination 
(recombination of an free electron with an ion). Grain charges should therefore be 
calculated in the network simultaneously. 

We consider neutral and singly charged dust grains. Gaseous species and dust grains 
interact through electron-grain, ion-grain collisions and adsorption and desorption of neutral 
species. Table 1 lists the generalized scheme of gas-grain interactions in our model. 

A collision between a negatively-charged grain and a positively- charged grain result 
in two neutral grains. The collision rate of a grain of radius a and charge le with another 
grain of radius a' and charge Ve is given by 



k 



l\2 



gg 



ir(a + a) 



8kT 
vr/i s 



1/2 



We 



'„2 



a')kT\ ' 



(2) 



where \x g is the reduced mass of the colliding grains ( lUmebayashi Ill983l ). We simply 
assume all grains have the same radius of 10~ 5 cm and density of 3 g cm -3 . In the case of 
electron-grain collisions, the rates are given by 



k 



tg 



S P ira z 



8kT 



Tvm f 



1/2 



J r 



akT 



-I 



(3) 



where m e is the mass of electron. An electron sticks to a grain upon collision with a 
sticking probability of S e ., wh ich strongly depends on the temperature. We calculate S e 
following lUmebayashi I (119831 ) . assuming a potential depth of the grain equals to 2 eV. J is 
a dimensionless reduced rate coefficient including the polarization effect of the grains by 
the electric field of app roaching species. We applied the approximate formulas given by 
Draine & Sutml Jl987f ). 

J(t,i/ = 0) = !+(-) , (4) 



J(t, v<0) 



V 






1 - - 






T- 







2u 



1/2 



(5) 



The collision of a cati on with a negatively charged grain is theoreti c ally d iscussed in 
Watson fc Salpeter I ( 119721 ) , and summarized by lUmebayashi fc Nakano I (ll980f ) . According 
to these references, there are two possible processes when a cation approaches a negatively 
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charged grain; a free electron tunnels to the cation before it hits the grain, or the 
recombination occurs on the grain surface. For simplicity, we assume the former case, in 
which the products and branching ratios of the recombination would be the same as the 
recombination in the gas phase. In the case of the collision of a cation with a neutral grain, 
recombination probably occurs on grain surfaces, and the products desorb to the gas phase, 
using the excess energy (T. Umebayashi 2012, private communication). Again, we assume 
that the products and branching ratios are the same as the corresponding recombination in 
the gas phase, for simplicity, although they might be different. For example, if X + + e~ — > 
Y + Z in gas-phase recombination, corresponding grain-surface recombinations are 

X+ + gr- -)• Y + Z + gr, (6) 
X + + gr -> Y + Z + gr + . (7) 

Our network includes recombination in the gas phase and on grain surfaces for all cations. 
Since we assume that each branching ratio for two fragment channels is only 5% (Section 
2.2.1), our approach tends to produce smaller molecules via the dissociation of large species. 
If we assume that, for example, the branching ratios for all product channels are equal or 
the collisions lead to simple charge exchange (ex. X + + gr~ — » X + gr), formation rates of 
large organic species are enhanced. In other words, our model gives the lower limit of the 
abundances of large organic species. We discuss how much the abundances of large organic 
molecules are enhanced if we assume the equally weighted branching ratios for grain surface 
recombination in Section 4.2. 

Recently, anions have been detected toward both the prestellar core, TMC-1 



fjBrunken et al. 1120071 ). and the protostellar core, LI 52 7 (ISakai et al. 1120071 ; lAgundez et al. 
20081 ). Our network includes H~, C~, 0~, S~, CN~, and OH - . Since interaction between 



anions and grains are unclear, we assume simple charge exchange, 

X" + gr+ -»• X + gr. (8) 
The rates of grain-ion (both cation and anion) collisions are expressed as 

a {8kT\ 1/2 ~ ( akT le\ 

where rrij and qj are the mass and the charge of ions, respectively. 

Gaseous neutral species c ollide with dust grains and stick to the grain surface. AW08 



and Ivan Weeren et al. I ( 120091 ) assumed sticking probabilities of 0.5 and unity, respectively, 
independent of the temperature. Since the sticking probability should be lower at higher 
temperatures, we assume that the sticking probability for neutral species is 

S n = [1 + 0.04(T + T d ) - 5 + 2 x 10~ 3 T + 8 x lO^T 2 ]" 1 , (10) 



- 10 - 



where T and are the gas and dust temperatures, respectively. We assume that T and Td 
are the s ame. Although Equation ffTUl) is originally derived for the sticking probability of a 
H atom ( Hollenbach fc Mckee 19791 ) by fitting the numerical calculations, we use it for all 
neutral species. For example, the value of S n is 0.83 at 10 K and 0.54 at 100 K. We might 
underestimate S n , because a loss of kinetic energy would be larger for collisions of heavier 
species with grains. 

After species are adsorbed onto grain surfaces, they can desorb to the gas phase again. 
We adopt the same desorption energies as iGarrod &: Herbst I (120061 ) , in which grain-surfaces 
are assumed to be covered with H2O ice. In addition to the thermal desorption, we consider 
two non-thermal desorption processes. One is the cosmic-ray induced desorption; grains are 
heated up to 70 K by the cosmic-r ay bombardment, and it al lows temporal evaporations of 
species with low binding energies (IHasegawa fc Herbst Ill993l ). The other is the chemical 
desorption; energy released by the exothermic grain-surface reactions is used to desorb the 
products. The efficiency of the desorption is parameterized by orkk, which is the ratio 
between the surface-molecule bond frequency and th e frequency at which the energy is lost 
to the surface. We assume orkk = 0.01 referring to IGarrod et al. I (120071 ); roughly 1% of 
the products are desorbed into the gas phase. 



2.2.3. Grain Surface Reactions 



Calculation of grain-surface chemistry is performed in a similar way to IGarrod fc Herbst 



(120061 ). The diffusion energy barrier is set to be a half of the desorption energy of each 
species. The diffusion rate is exponentially dependent on the diffusion barrier divided 
by the grain temperature, and surface r eaction rates are determ ined by the sum of the 
diffusion rates for two reaction partners (IHasegawa et al. Ill992l ). If surface reactions have 
activation energy barriers, the barriers are overco me thermally, or via quantum tunneling, 



whichever is faster. The modified rate method (jCaselli et al. I Il998l ) is used only for 



reactions involving atomic hydrogen. Our model also includes photodissociation on grain 
surfaces by interstellar and cosmic-ray induced UV photons. 

Our model assumes that molecules are physorbed, which is appropriate as long as 
grain surfaces are covered with ice. At T > 150 K, however, ice mantle is subli mated, 
and chemisorption might become important (e.g., Section 4.2.2 of iTielens II2005I ). Since 
grain-surface reactions via chemisorption is out of the scope of this paper, we switch off 
adsorption of neutral gas-phase species onto grain surfaces, desorption of ice-mantle species 
to the gas phase, and grain surface reactions, if the temperature becomes greater than 200 
K. 
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2.2.4- Initial Molecular Abundances and Ionization Rates 



Table 2 lists the initial abundances in our model. The elemental abundances adopted 
here are called "lo w metal" values beca use of their strong depletions of species heavier 
than oxygen (e.g., 



Graedel et al. Ill982[ ). Species are initially assumed to be atoms or 



atomic ions except for hydrogen, which is in molecular form. We integrate the initial 
abundances for 10 5 yr with the typical conditions of dense molecular cores, T = 10 K, 
% = 1 x 10 4 cm -3 and visual extinction of Ay = 10 mag, which determine the abundances 
of our initial Bonnor-Ebert like sphere. Before the collapse begins, we further integrate 
the updated abundances under the initial hydrostatic conditions during 1.6 x 10 5 yr (i.e., 
sound crossing time of our core model), in order to obtain the starting abundances for the 
collapse, implicitly assuming that the core keeps its initial hydrostatic structure supported 
by turbulence. Our approach is similar to AW08 and van Weeren et al. (2009). 



In the hydrodynamic simulation, we calculate column densities of gas in the x—, y— 
and z— directions, E x , E y and E z , from the core outer edge to a position of each fluid 
parcel. We choose the minimum of E^, E y , E z and convert it to Ay, via the formula 
A v = (/ H E///m H ) x (5.34 x 10~ 22 cm 2 ) mag faohlin et al.~lll978l : ICardelli et al~lll989h . 



where /h is the mean number of hydrogen nuclei per species (1.67), fi is the mean molecular 
weight (2.3) and mn is the mass of a hydrogen atom. For the fluid parcels we follow, the 
initial values of Ay are greater than 10 mag, and Ay increases as the collapse proceeds. 
Hence the interstellar UV radiation has little effect on chemistry in our model. 

Cosmic rays are the main source of ionization in molecular clouds. The cosmic-ray 
ionization rat e of hydrogen is set to be 1.3 x 10~ 17 s -1 . The attenuation length of cosmic-ray 



is 96 g cm -2 ( lUmebayashi fc Nakano Ill98ll ). which is much larger than the column density 
of typical molecular clouds, but can be smaller than the column density, E, calculated 
above. The decay of short-lived radio active isotopes, wh ich releases energetic particles , 
also ionize the ISM with rates f = (7.6-11) x 10~ 19 s" 1 dlJmebavashi k Nakanoll2009l ). 
Thus for direct cosmic-ray ionization and cosmic-ray-induced photo reactions, we assume 
that the effective cosmic-ray ionization rate is expressed as 



i eff = 1.3 x 10~ 17 exp 



96 gcnr 



1 x 10 



-18 



'ID 



In our model, the decay of radio active isotopes is the main source of ionization within 
the first core, while the cosmic-ray dominates at the surface of the first core and infalling 
envelopes. 
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3. RESULT 

3.1. Physical Structure of the Core 

Figure [3] shows the spatial distributions of the gas number density, radial velocity (v r ), 
and temperature perpendicular to (a-c) and along the rotational axis (d-f). We averaged 
the distributions perpendicular to the rotational axis in the z = plane in the azimuth 
direction. The solid lines represent the distributions at £fc — 1500 yr, while the dashed 
lines represent the distributions at tpc — 2800 yr. We define r and R as \J x 2 + y 2 + z 2 
and \/ x 2 + y 2 , respectively. Although the core has the rotation, the core structures are 
almost spherical at r > 10 AU both in the middle and late stages. The density has a 
power law distribution of r~ 2 at r > 100 AU and r~ L5 at 10 AU < r < 100 AU. The 
effect of the rotation becomes apparent in the inner region (r < 10 AU) especially in the 
late stage, in which infalling gases have larger angular momentum. The accretion onto the 
first core is supersonic, and the accretion shock occurs at the first core surface. The first 
core is quasi-hydrostatic, and in the middle stage its radius is ~6 AU perpendicular to the 
rotational axis and ~5 AU along the rotational axis, respectively, while it is ~10 AU and 
~3 AU in the late stage. 

The temperatures both in the envelope and the first core increase with time as the 
first core grows, except for the low-temperature region in 5 AU < R < 10 AU around the 
equatorial plane. This region experiences strong expansion perpendicular to the rotational 
axis in a dynamical timescale due to the angular momentum transport. Cooling by the 
expansion is more efficient than the radiation heating, since the gas is optically thick. 

3.2. Spatial Distribution of Fluid Parcels 

Figure 0] shows the spatial distributions of fluid parcels in the middle (left) and late 
(right) stages. Each parcel is color-coded by its initial distance from the core center (above 
the equatorial plane) or its initial specific angular momentum (below the equatorial plane). 
In the middle stage, the first core is composed mainly of the fluid parcels which were initially 
at r ~ 1150-1450 AU. In the late stage, these parcels have moved to the midplane of inner 
radii due to contraction, and the outer layers of the first core is composed of parcels falling 
from the initial radii of r ~ 1700-1850 AU. It indicates that the accretion of the envelope 
onto the first core is active, and the gas inside the first core is continuously compress into 
the center. Parcels closer to the rotational axis have smaller initial angular momentum. 
While parcels fall almost radially in the envelope, they move in the z direction in the first 
core rather than radially until they reach the equatorial plane. Near the equatorial plane, 
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parcels move inward in the direction perpendicular to the rotational axis by losing the 
angular momentum. 



3.3. Molecular Evolution in a Single Fluid Parcel 

Here, we show the thermal and chemical histories of a single fluid parcel, which are 
helpful to understand the spatial distributions of molecules presented in Section 3.4. Figure 
[5] shows the trajectory of the parcel, which is initially at r = 1300 AU and 70° degrees from 
the rotational axis. The first core surface (white dashed line) is defined as the isoline of the 
minimum gas density in the region where the local sound speed is greater than the radial 
velocity. Figure [6]^a) shows the temporal variations of the position of the parcel in r, R, z. 
The horizontal axes is set to be t fina i — t, where tfi na i represents the final stage of our model; 
£fc = corresponds tfi na i — t ~ 3000 yr. The vertical dotted lines in Figure [6] represent the 
time when the parcel passes the accretion shock at the first core surface (tfmai —t— 1930 
yr). The parcel initially infalls almost spherically. The effect of the rotation becomes 
apparent inside the first core, and the trajectory becomes helical. The parcel falls to the 
equatorial plane within 10 2 yrs after entering the first core, while it moves inward in the 
direction perpendicular to the rotational axis over 10 3 yrs by losing its angular momentum. 

Figure MJo) shows the temporal variations of the temperature in the parcel. We classify 
the history of the parcel into three phases in terms of the temperature; cold phase (10-20 
K), warm-up phase (20-200 K) and hot phase (>200 K). Dominant chemical processes 
vary among the three phases. In the cold phase, adsorption onto and hydrogenation 
on grain surfaces are dominant, while in the warm-up phase, the molecular evolution is 
caused by grain-surface reactions among species with heavy elements and evaporation 
of ice mantle species. In the hot phase, the molecular evolution is caused by gas phase 
reactions. This classification has been used in the chemistry in protostellar stage (e.g., 



Herbst fc van Dishoeck 1120091 ). It should be noted, however, that in the first core stage 
the density is higher and the duration of the warm-up and hot phases are shorter than in 
protostellar stage. Table 3 summarizes the physical conditions and duration of the three 
phases in the parcel. 

In the following, we show the molecular evolution in these three phases. However, the 
temporal variations at 20 K < T < 1000 K is difficult to see because the temperature rises in 
a short time scale (~10 3 yrs). In order to zoom-up the evolution at that temperature range, 
each panel in Figures |6]^c-h) is divided into three phases based on the above classification. 
The panels (c), (d) and (e-h) show the number density of H nuclei, temperature, and 
molecular abundances relative to H nuclei in the parcel, respectively. The physical 
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conditions as well as the molecular abundances look discontinuous at the boundaries of the 
phases, which is due to the different scaling of the horizontal axis. We assume that the 
parcels stay for 1.6 x 10 5 yr in its initial position, as mentioned in Section 2.2.4, but this 
pre-collapse phase is omitted in Figure El 

3.3.1. Cold Phase 

When the core begins to contract (ifi na i — t = 7 x 10 4 yr), molecules with heavy elements 
are mostly exists on grain surfaces. The main carbon reservoirs are CO, CH4 and H 2 CO in 
order of high to low abundances. Water, which is formed by hydrogenation of oxygen atoms 
on grain surfaces, is as abundant as ~10~ 4 . The nitrogen mainly exists as N2 and NH3. All 
these molecules mostly exist on grain surfaces. 

At the beginning of the contraction, the parcel migrates inward almost spherically 
and isothermally. As the parcel falls to inner denser regions, the abundances of gaseous 
molecules decrease due to freeze-out. Hydrogenation on grain surfaces occur subsequently, 
and the abundances of formaldehyde (H2CO) and methanol (CH3OH) slightly increase. At 
the end of the cold phase, the abundances of CH 3 OH and H 2 CO are 3 x 10~ 6 and 1 x 10~ 5 , 
respectively. 

3.3.2. Warm-up Phase 

After the first core is formed, the parcel continues to infall, and the temperature in 
the parcel starts to rise. The parcel passes the accretion shock at the first core surface 
at ^finai — t = 1930 yr, when the number density of H nuclei is 3 x 10 10 cm -3 and the 
temperature is 80 K (vertical dotted lines in Figure E]). Since the shock heating layer is not 
spatially resolved in our hydrodynamic simulation, our chemical model neglect this shock 
heating. We discuss the effect of the accretion shock on chemistry in Section 4.4. 

During the warm-up phase, species are desorbed into the gas phase when the 
temperature reaches their sublimation temperatures. CO is desorbed to the gas phase when 
the temperature reaches 25 K, which is slightly higher than the typical interstellar value 
of 20 K, because of the higher density (~10 9 cm -3 ) of the envelope in the first core stage. 
On the other hand, species with the higher binding energies are not desorbed yet; they 
migrate on grain surfaces and react with each other at the lukewarm temperature of T ~ 
several 10 K. AW08 found that large organic molecules, such as HCOOCH 3 and CH3OCH3 
are efficiently formed at T = 20-40 K, and reach the terminal abundances of ~10~ 9 in 
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the protostellar core. In our model, HCOOCH3 is formed at T = 30-50 K, but reaches 
abundances of only 2 x 10~ n . The abundances of formic acid (6 x 10 -11 ) and CH 3 CN 
(6 x 10 -9 ) change little in the warm-up phase. The formation of large organic species in 
the warm-up phase is inefficient in our model because of the short time scale spent in the 
lukewarm region, which is determined by the size of the lukewarm region divided by the 
infalling yelocity: it is only ^10^ yrs in our model, while it continues ~10 4 yrs in AW08 
(see also iGarrod &: Herbst II2006I ). These large organic molecules start to be desorbed at 
T > 100 K. 



3.3.3. Hot Phase 

In the hot phase, ice mantles of grains are completely evaporated, and gas phase 
reactions determine the molecular abundances. The temperature and density of H nuclei 
reach 2000 K and 8 x 10 15 cm -3 , respectively at the final step. The duration of the hot 
phase is ~2 x 10 3 yrs. 

At the temperature of 200-500 K, the composition stays almost the same as in the 
warm-up phase. At T > 500 K, most molecules are destroyed, and simple molecules are 
reformed. For example, H 2 CO starts to be destroyed by collisional dissociation when 
the temperature reaches ~500 K. CH 3 OH is also collisionally dissociate at T ~ 1200 K. 
Methane, which is the second carbon reservoir, is converted to CO by neutral-neutral 
reactions in several hundred yrs. As a result, the abundance of CO increases by a factor of 
2. The abundance of H 2 0, which is the main oxygen reservoir, remains the almost constant 
value of ~10~ 4 , which is mostly determined in the precollapse phase. The dominant 
reservoirs of the nitrogen are N 2 and NH 3 . At T > 1600 K, NH 3 is converted to N 2 through 
the reactions of NH 3 + H ->■ NH 2 + H 2 , and NH 2 + NO N 2 + H 2 on the time scale 
of ~10 yrs. The former reaction has the potential energy barrier of 5.2 x 10 3 K. The 
abundance of hydrogen atom increases with time mainly by collisional dissociation of H 2 . 
When the temperature reaches ~2000 K, ~0.1% of hydrogen nuclei is in hydrogen atom. 
As the abundance of hydrogen atom increases, OH is formed through the reaction, H 2 + 
H — y OH + H 2 , which has the potential energy barrier of 9.7 x 10 3 K. 

Some larger molecules are temporary formed in the gas phase. Methyl formate and 
CH 3 OCH 3 are formed by the grain-surface recombination of HCOOCHj and CH 3 OCHj~, 
respectively, despite the small branching ratio of 5%. These formation paths are efficient 
until CH 3 OH, which is the precursor of HCOOCH^ and CH 3 OCHj, is destroyed by 
collisional dissociation. The peak abundances of HCOOCH 3 and CH 3 OCH 3 are 5 x 10~ n 
and 9 x 10~ n , respectively. At the temperature of ~700 K, acetylene (C 2 H 2 ) is formed 
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by the reaction of H2 and C3H4, which is formed on grain surfaces during the cold phase. 
The reaction has the potential barrier of 2.3 x 10 4 K. The abundance of C2H2 reaches 
~10 -6 at the maximum. C2H 2 is destroyed through the reactions of C2H2 + OH — > CH 2 CO 
+ H, and CH 2 CO + M — > CH 2 + CO, where M represents the third body (mainly H 2 ). 
These reactions have the potential energy barrier of 5.0 x 10 3 K and 3.9 x 10 4 K. Other 
unsaturated carbon-chain molecules, C x H y , are temporally formed via the neutral-neutral 
reactions, C^H^! + H — > C x H y + H 2 , which have the potential energy barrier of several 10 3 
K. HCN is converted to HC 3 N through the reactions of HCN + H CN + H 2 , and C 2 H 2 
+ CN — > HC3N + H. The former reaction has the potential energy barrier of 1.2 x 10 4 K. 
Larger cyanopolyynes (HC n N, n — 5,7, 9) are also formed by the subsequent neutral-neutral 
reactions of HC 3 N + C 2 H, HC 5 N + C 2 H and HC 7 N + C 2 H. The terminal abundances of 
HC„N (n = 3, 5, 7, 9) are 4 x 10" 7 , 1 x 1(T 7 , 2 x 1(T 8 and 4 x 10" 9 , respectively. These 
cyanopolyynes survive until the end of the simulation, because the neutral-neutral reactions, 
C n N + H 2 — > HC n N + H, which have the potential barriers of 2000 K, are compensated 
by the formation via the collisional dissociation reactions, HC n N + M — > C n N + H + M. 
Therefore, the abundances of cyanopolyynes depend on the products of the collisional 
dissociation reactions assumed in our model. 

In summary, the total molecular abundances (gas and ice combined) of many species 
in the cold phase mostly remain unaltered until the temperature reaches 500 K. Above 500 
K, the molecules with high abundances start to be destroyed, and are reformed to simple 
molecules, such as CO, H 2 and N 2 . Some molecules are formed in the warm-up phase 
and the hot phase; HCOOCH3 is formed on grain surfaces during the warm-up phase, 
while it is formed in the gas phase during the hot phase. Carbon-chains, such as C2H2 and 
cyanopolyynes, are effectively formed in the gas phase in the hot phase. 

3.4. Spatial Distribution of Molecular Abundances 

In order to derive the spatial distributions of molecular abundances, we calculate 
molecular evolution in 10 parcels, which are selected from the 10 5 parcels traced in our 
RHD simulations. Figure [7] shows the distributions of the parcels, which are color-coded 
according to the fractional abundances of gaseous CH3OH. The contour lines and gray 
scales depict the temperature and number density, respectively, in logarithm scale. The 
abundance of CH 3 OH jumps at r ~ 5 AU in the middle stage and r < 10 AU in the late 
stage, where the temperature reaches its sublimation temperature (~100 K). At the outer 
radii, CH3OH exists on grain surfaces. The spatial distributions of other many molecules 
are similar, while sublimation radii vary among species. In the late stage, the abundance of 
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CH3OH drops in the central region (r < 1 AU), because CH3OH is collisionally dissociated. 

Figure [7] indicates that the distributions of CH 3 OH abundance follows the local 
temperature. To confirm this, we plot the abundance of CH 3 OH as a function of 
temperature in the middle stage (gray square) and the late stage (black cross) in Figure [HJ 
It clearly shows that the CH3OH abundances follow the local temperatures. Note, however, 
the plots in the middle stage and the late stage are not identical. It indicates that the fluid 
parcels with the similar temperature at each evolutionary stage have experienced similar 
thermal histories and molecular evolutions. At the temperature of around 100 K and 1200 
K, the CH3OH abundances are scattered. The CH3OH abundances are determined by 
the sublimation at ~100 K and collisional dissociation at ~1200 K. The rates of these 
processes have the exponential dependence on the temperature, and a little difference of the 
thermal histories may make the scattering. In the late stage, there are two branches at the 
temperature below ~100 K. The lower branch represents the CH3OH abundances in the 
local low-temperature region inside the first core, which is mentioned in Section 3.1, while 
the upper branch represents the abundances in the envelope. Since the density is much 
higher in the first core than in the envelope, the degree of the depletion is much sever in the 
first core. 

The spatial positions of the parcels at each snapshot constitute an irregular grid for 
which the molecular abundances, the density and the temperature are known. We have 
interpolated the molecular abundances in the density-temperature space with the nearest 
neighbors method, and obtained the table of the molecular abundances as the function 
of the density and the temperature at each snap shot. These tables are used to derive 
the molecular abundances at each evolutionary stage at the points of three-dimensional 
regular Cartesian-grids for which the temperature and the density are known from our 
hydrodynamic simulation. We use these data to derive the profiles of molecular abundances 
and column densities. Since the molecular evolution in the star-forming regions is a 
non-equilibrium process, in general, the abundances are not given as a function of only the 
local physical conditions, but also depend on the past physical conditions. But in our case, 
the abundances can be given as a function of the local temperature and density at each 
evolutionary stage, because the flow pattern in our hydrodynamic simulation is simple, 
and again the fluid parcels with the similar physical conditions at each evolutionary stage 
have experienced the similar thermal histories and molecular evolutions. 

Figure [9] shows the radial distributions of the physical parameters (a, b) and molecular 
abundances (c-h) in the z = plane in the middle stage (a, c, e, g) and the late stage 
(b, d, f, h). The distributions along the rotational axis are shown in Figure [TOl which is 
available in the online version of this paper. Since the temperature distributions are almost 
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spherical (Figures EJjc) and[2](f)), and since the molecular abundances depend mainly on 
temperature, the molecular distributions are similar in Figure |9] and Figure [101 except near 
the first core (r < 10 AU). In the remainder of this subsection, we discuss the abundance 
profiles in detail. 



3.4- 1- Simple Molecules 

Simple molecules are often observed both in prestellar and protostellar cores, and are 
used to probe the physical conditions. Among the neutral species observable in radio, CO 
has the lowest binding energy, and thus its sublimation radius is relatively large; r ~ 40 
AU in the middle stage and r ~ 100 AU in the late stage. It could be observed as a CO 
emission spot inside the heavy CO depletion zone. Sublimation radii of U2CO and HCN 
are also relatively large; r ~ 15 AU in the middle stage and r ~ 30 AU in the late stage. In 
the late stage, there is the low-temperature region (T ~ 10 K) inside the first core on the 
equatorial plane (Section 3.1). All gaseous species including CO are on grain surfaces there, 
because of the low temperature and high density. The freeze-out time scale is estimated 
to be only hours. Sublimation radii in the equatorial plane of the polar species, such as 
H 2 and NH 3 are small: R ~ 5 AU. Since the temperature distribution deviates from the 
spherical system in the small scale (r < 10 AU), the sublimation radii of the polar species 
are slightly larger along the rotational axis; z ~ 8 AU. In the middle stage the abundance 
of CH 4 is constant, but in the late stage the abundance drops at r ~ 1 AU. The peak 
abundances of species mentioned above do not change much with time. 

The dominant positive ion is H^" outside the sublimation radius of CO. At the inner 
radii, the dominant positive ion varies in accordance with proton affinity; HCO + dominates 
inside the CO sublimation radius, but it is replaced by H 3 CO + and then NH4 , as H 2 CO 
and NH3 sublimate subsequently. The ionization degree drops toward the first core, due to 
increasing density. Dominant charge carrier is the charged grains at n > 10 10 cm -3 . 



3.4-2. Carbon-Chain Molecules 



Unsaturated carbon-chain molecules are typically associated with cold dense cloud 



cores, but have also be en detected towar d the Class-0/I sourc es, L1527 fjSakai et al. 1120081 ) 



and IRAS 15398-3359 flSakai et al. 1120091 ) . ISakai et al. I (120081 ) suggested that carbon chains 



in the protostellar cores are formed from the sublimated CH 4 in the lukewarm region near 
the protostar. CH 4 reacts with carbon ion to form C2H^", which is a precursor of larger 
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unsaturated carbon-chains. This formation scenario of carbo n chains in the prot ostellar 



cores was named warm carbon-chain c hemis try (WCCC) by ISakai et al. I (120081 ) . and 



confirmed by AW08 and lHassel et al. I ( 120081 ) . 



We found WCCC does not efficiently occur in the first core stage; unsaturated 
carbon-chains do not increase at sublimation radii of CH4, r ~ 90 AU and 30 AU in the 
middle and late stages. The abundance of carbon ion is lower due to the higher density 
of the envelope than in protostellar stages. Instead high temperature chemistry produces 
C2H2 and cyanopolyynes at r < several AU. In the late stage, the abundance of C2H2 drops 
at r < 1 AU, where longer cyanopolyynes, HC n N (n = 5, 7, 9), increase. 



3.4-3. Large Organic Molecules 

Large organic molecules sublimate at several AU. The total of gas and ice abundances 
of each species are almost constant throughout the core and envelope except for HCOOCH3 
and CH3OCH3. HCOOCH3 is partly formed on grain surfaces at the lukewarm region, 
which expands from several tens AU (middle stage) to ~ 100 AU (late stage). CH3OCH3 
is mainly formed in the gas phase at r < 2 AU. In the central region, the abundances of 
CH3OH, HCOOCH3 and CH3OCH3 drop in the late stage. 



4. DISCUSSION 

4.1. Molecular Column Densities 

Detection and observation of first cores is one of the most important challenges in 
ALMA era. First cores are characterized by a AU-size hydrostatic core and low-velocity 
outflow with a wide opening angle. While size of the core can be measured by high angular 
resolution dust continuum observation, molecular line observations are also necessary to 
confirm the kinetic structure. Molecular abundances obtained in our model are useful to 
predict where the various molecular lines would arise. We especially focus on large organic 
species, which are abundant at th e vicinity and the surface o f the first core. Since detailed 



radiation transfer modeling (e.g., iTomisaka fc Tomida 1 1201 if ) is out of the scope of this 



work, here we compare the column densities of assorted species in our model with those 
estimated from the observations towards a prototypical hot corino, IRAS 16293-2422, to 
evaluate if our model column densities of large organic species are high enough to produce 
lines. 
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Figure [TT] (solid black line) shows the cumulative column densities of gases along 
the rotational axis from the outer edge of the core to a height z above the midplane, 
N z (z) = J'^ndz'. The vertical solid lines represent the first core surface. The vertical 
dashed and dotted lines represent the heights where the optical depths of dust continuum 
reach unity at the wavelength of 100 /xm and 1 mm, respectively. The optical depth of dust 
continuum, r, is estimated as 



T 



8 x 10 25 cnr 



A 



1 mm 



(12) 



referring to Hildebrand (1983). At 100 jum the optical depth reaches unity in the envelope. 
At 1 mm, on the other hand, we can look into the envelope and the surface of the first 
core, although the optica l depth exceeds unity at deeper layers of the first core (see also 
Saigo fc Tomisaka Il201l[ ). The gases in front of the r lmm ~ 1 plane can be probed by 
absorption lines, since the gas in front is cooler than the Ti mm ~ 1 plane. In the line of sight 
offset from the rotational axis where dust continuum is optically thin (R>6 AU), the line 
could be emission (Figure [T2l) . We calculate the molecular column density in front of the 
r imm = 1 plane in the face-on system, 



iV, 



rii(x,y,z')dz', 



(13) 



z T =i{x,y) 



where z T= \ represents the height z where the optical depth reaches unity at 1 mm, and n 
represents the number density of species i. Figure [T2l shows N Zr=l i of assorted molecules £ 
a function of R in the middle stage (upper panel) and the late stage (lower panel). 



At R < 6 AU the first core is optically thick for 1 mm dust continuum, and N Zi _ =lti is 
almost constant. N Zt=1: qo sharply increases inwards at its sublimation radius (several tens 
AU) and at the first core surface (several AU). N Zr=lt i of large organic molecules sharply 
increases at their sublimation radii which is close to the first core surface. The peak value 
of N Zt=1 ^ for each species are not much different between the middle and late stages. 
CH3OCH3 is an exception; its peak value of A' 2t=1! ch30CH3 increases with time. 

The emission of thes e complex organic molecules have been detected towards severa l 



low mass protostars (e.g., ICazaux et al. 1120031 ; iBottinelli et al. 1 12004 ; iMaret et al. II2005I ). 



IRAS 16293-2422 are A^ C h 3 oh = 1.1 X 10 1 8 cm" 2 , Nncnn™ , 



2004h . and AT C h 3 och 3 = 6.4 x 10 15 cm" 2 flChandler et al. 



x 10 15 cm" 2 ( 


Kuan et al. 


Kuan et al. 


(2004 


) also 



derived Ai3ch 3 oh = 8.1 x 10 cm" " 2 . It indicates that the 12 CH 3 OH emission is optically 
thick, since the carbon isotope ratio of CH3OH, Ni2c,k 3 cm / Ni3r,K 3 nn ~ 14, is much smaller 



than the elemental abundance ratio, [ C/ C] ~ 60, in ISM (e.g.. lLanger fc Penzias 111993 ; 
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Lucas fc Liszt II1998I ). Our model column densities of CH3OH and HCOOCH3 in the middle 
and late stages at R < 6 AU are comparable or larger than those values. Hence it would 
be reasonable to expect that the column densities of CH 3 OH (and HCOOCH3) at R < 6 
AU in our model are high enough to be 77 ~ 1 at least in some transitions. Since the gas 
in the first core surface and the envelope is cooler than the first core, we expect absorption 
lines, if we observe the object from the pole-on angle. Assuming the lines and continuum 
are thermalized, the expected difference of brightness is 

AT b = \T bg - T fg \[l - exp(-r,)]> (14) 

where T bg and Tf g are the temperatures at Ti mm = 1 and 77 = 1, respectively. We assume 
that 77 for CH3OH reaches unity at the height where the model column density of CH3OH 
reaches 1.1 x 10 18 cm" 2 . In the middle stage for CH 3 OH, AT b is >40 K at R < 6 AU, 
which is detectable in 1.6 hours toward the distance of 120 pc with S/N > 10 and velocity 
resolution of 0.1 km/s in the extended configuration (~0.1") of ALMA at Band 7 (345 
GHz). More integration time would be required to detect HCOOCH3, because the 77 = 1 
plane for HCOOCH3 is closer to the ri mm = 1 plane than that for CH3OH is (i.e., AT?, for 
HCOOCH3 is smaller) (see Figure [XT]) . We found that the column densities of HCOOH 
and CH 3 CN in our model are also comparable to the observed values in the hot corinos. 
Although, detailed radiation transfer calculation is needed to confirm these estimates. 

On the other hand, the column density of CH3OCH3 in our model is much less than 
the observed toward IRAS 16293-2422 both in the middle and late stages, and thus would 
not be high enough to produce observable lines. 

Our results show that first cores and the surrounding warm envelopes could be observed 
as very compact hot corinos without stellar signatures, and the large organic molecules can 
be good tracers of first cores. In the rest of this paper, we call the warm gases with the 
large organic molecules a prestellar hot corino. 



4.2. Uncertainties in the Reaction Network Model 

There are several caveats to our model abundances and column densities of large 
organic molecules presented above. As mentioned in Section 3.3.3, H 2 CO is collisionally 
dissociated at lower temperatures (~500 K) than other large organic molecules, since the 
value of 7 (1.75 x 10 4 K) is small for the collisional dissociation, H 2 CO + M — > H 2 + CO 
+ M, compared with those of other large organic molecules. It should be noted that the 
value of 7 in the collisional dissociations listed in Table A2 are basically determined by the 
difference between the formation enthalpy of reactants and products, but the products and 
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branching ratios for various possible sets of products are not well constrained by laboratory 
experiments (Section 2.2.1). If there exists the channels with lower 7 for the destruction 
of other large organic molecules, the abundances of these molecules at T > 500 K could 
decrease by several orders of magnitude relative to our results. The column densities of the 
large organic molecules in the middle stage (Figure [121 upper panel), however, would not 
change, since the temperature is lower than 500 K at Ti mm < 1. The column densities in 
the late stage (Figure [121 lower panel) would decrease by about one order of magnitude. 

The abundances of HCOOCH3 and CH3OCH3 in our model depend on the branching 
ratios of grain surface recombination, which contributes to their formation in the hot phase. 
We recalculate the molecular abundances in the representative fluid parcel discussed in 
Section 3.3, assuming all the product channels of grain surface recombination (both with a 
negatively charged grain and with a neutral grain) are equally weighted. Then, the peak 
abundances of HCOOCH3 and CH3OCH3 are enhanced by a factor of 8 and 12, respectively. 
The column densities would be enhanced almost by the same factor, since grain surface 
recombination is effective around/in the dense first core. In contrast, the abundance of 
CH3OH does not strongly depend on the assumed branching ratio of the recombination, 
since CH 3 OH is mostly formed by hydrogenation on grain surfaces in the cold phase. 

While the cosmic ionization rate is set to be 1. 3 x 10~ 17 s~ l in our model, a larger 



value of 5 x 10~ 17 s" 1 is suggested by recent work f IDalgarno II2006I ). We recalculate the 
molecular abundances in the representative fluid parcel with the cosmic-ray ionization rate 
of 5 x 10~ 17 s^ 1 . The peak abundances of large organic molecules are enhanced by a factor 
of 2-5 compared to our original model with the cosmic-ray ionization rate of 1.3 x 10~ 17 s" 1 . 

In summary, our model might underestimate the abundances of large organic molecules 
at T < 500 K by one order of magnitude, while our model might overestimate them 
at T > 500 K by several orders of magnitude. The model column densities might be 
underestimated in the middle stage by one order of magnitude, while those in the late stage 
could be higher or lower by one order of magnitude. 



4.3. Comparison with Protostellar Hot Corinos 



Here, we compare the physical structures and chemical compositions of the prestellar 
hot corino (prestellar HC) with those of a protostellar hot corino (protostellar HC). W e refer 
to the spherical RHD model of a star- forming core by iMasunaga fc Inutsuka I (|2000[ ). and 
AW08 for the protostellar hot corino model. Such comparison is useful to unveil differences 
between prestellar and protostellar HCs. It should be noted, however, that the physical 
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structures of the cores are different; e.g., while the models of iMasunaga Sz Inutsuka I (j2000f ) 



and AW08 are non-rotating cores of 3.9 Mq, our model is the rotating core of 1 Mq. 

The top two columns of Table 4 show the mean radius and the gas number density 
of T = 100 K region. It is clear that the prestellar HC is more compact and denser than 
the protostellar HC by an order of magnitude; the protostar is more luminous and their 
envelopes are more sparse than in the first core stage. 

We compare the abundances of the large organic molecules in the central region 
(r ~ 5 AU) of our model with th ose of the protostellar HC mo del. While AW08 traced 13 



parcels in the RHD simulation of IMasunaga fc Inutsuka I (l2000f ) to calculate the molecular 



abundances, the reaction network used in AW08 is different from our network. In order to 
make comparisons easy, we revisit the calculation of AW08 with some modifications; we 
calculate the temporal variations of the molecular abundances using our network model in 
the parcel which reaches r = 15 AU at the protostellar age of 9.3 x 10 4 yr. We adopt this 
parcel, because the abundances are mostly constant at r < 100 AU according to AW08. 
While AW08 used the initial temperature of 6 K, we set it at 10 K when the temperature 
is lower than 10 K in the original model. Table 4 shows the abundances of the large 
organic molecules at 15 AU in the protostellar HC model and 5 AU in the prestellar HC 
model. Comparing the second row of Table 4 with Table 3 of AW08, our modifications 
in the network and temperature lead to order-of-magnitude changes in the abundances of 
listed species except for HCOOCH 3 . Among the listed species, the abundances of H 2 CO 
and CH3OH are similar in the prestellar HC model and the protostellar HC model. These 
species are formed during the cold phase by the subsequent hydrogenation reactions of CO 
on grain surfaces. The other species are less abundant in the prestellar HC model. In the 
protostellar HC model, they are efficiently formed in the lukewarm region. On the other 
hand, in the prestellar HC model the lukewarm region is small, and the time to form these 
species is limited. It indicates that the abundances of these molecules are sensitive to the 
luminosity of the central object. 



4.4. Accretion Shock Heating Layers 



The supersonic accretion of the envelope material to the first core causes the shock 
wave. This accretion shock is isothermal, because the radiation from post-shock gas heats 
up pr e-shock gas, and the temperature in the pre-shock and post-shock gases become equal 
(e.g., lOmukai 1 120071 : ISaigo fc Tomisaka Il201ll ; ICommercon et al. Il201ll ). Between the pre- 
and post-shock gases, the shock heating layer is formed, in which the gas temperature 
reaches ~2000 K at the maximum. In contrast, the grain temperature increases to only 



-24- 



200 K. The thickness of this layer is ~10 3 AU, and the cooling time scale of the gas 



is ~1 hour after crossing the shock front (see Fig. 3 of ISaigo fc Tomisaka II2011I ). In our 
hydrodynamic simulation the shock heating layer is not spatially resolved, and thus our 
chemical calculations neglect the temperature rise in the shock heating layer. In order to 
evaluate the effect of shock heating on chemistry, we recalculate the molecular abundances in 
the parcel discussed in Section 3.3 by adding 1 hour integration at the constant temperature 
when the parcel reaches the first core surface. Since the temperature of the shock heating 
layer should vary depending on the mass of the first core and the radius where the shock 
occurs, we calculate two cases: the gas temperature of 1000 K and 2000 K. The density is 
assumed to be the same as the post-shock gas along the trajectory, riu = 8.5 x 10 11 cm -3 . 
The grain temperature is assumed to be 200 K, at which temperature ice-mantle species are 
almost evaporated. In the pre-shock and post-shock regions, on the other hand, the species 
with the high binding energies, such as CH 3 OH, is mostly freeze-out on grain surfaces, since 
the grain temperature is about 80 K. To focus on high temperature chemical processes, we 
discuss changes of the total (gas and ice combined) abundances in this subsection. 

In the case of the gas temperature of 1000 K, the shock heating does not affect the 
chemistry. Even in the case of the gas temperature of 2000 K, the heating has a limited 
effect on the chemistry; among abundant species (nj/nn > 10~ 12 ), only 23 species increase 
or decrease by a factor of >2 in the shock heating layer. Figure [13] shows the spatial 
variation of the temperature (upper panel) and molecular abundances (lower panel) along 
the trajectory after the parcel experience the shock heating of 2000 K at r ~ 7 AU and 
^finai — t= 1930 yr (see Figure E]). The horizontal axis in Figure [T31 is the position of the 
parcel in r. In the lower panel, the solid lines show the molecular abundances with shock 
heating, while the dashed lines show our original model without the shock heating. The 
abundances of HCN and CH 3 OH are almost the same between the two models. H 2 CO is the 
most sensitive to the shock heating; it decreases approximately five orders of magnitude by 
collisional dissociation to form CO and H 2 . Since the abundances of CO and H 2 are higher 
than that of H2CO, the effect of H2CO dissociation to other molecules is insignificant. 
The difference between the two models diminishes when the temperature reaches ~500 
K, at which collisional dissociation of H 2 CO becomes efficient in the original model. 
Carbon species, such as C 2 H 2 , are also sensitive to the shock heating; their abundances are 
temporally enhanced by approximately four orders of magnitude by the reaction of C3H4 
and H 2 with the potential barrier of 2.3 x 10 4 K. The difference in the C 2 H 2 abundance 
diminishes when the temperature reaches ~800 K. Therefore our results in Section 3 might 
overestimate the abundance of H 2 CO, while underestimate the abundance of C 2 H 2 in the 
first core except for the hot central region (r < 3 AU). 

While our simple calculation shows that the shock heating affects abundances of some 
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molecules, the effect is limited and the shock heating does not lead to a significant change 
in our results presented in Section 3. 



4.5. Toward the Chemistry in Circumstellar Disks 

The first core is a transient object formed prior to the protostar. Here, we briefly 
discuss what happens after the formation of the protostar and the circumstellar disk in 
terms of chemistry. 

There are several studies which investigate the formation and evolution of the 
circumstellar di s k using thre e-dimensional hydrodynamic simulations. We refer to 



Machida et al. I (120101 120111 ). which investigated a longest-term evolution of the 
circumstellar disk, starting from the gravitational collapse of the molecular cloud core. 
While the central part (r < 1 AU) of the first core collapses to form a protostar, the outer 
part directly evolves to the cicumstellar disk with a Keplerian rotation. In that case, 
the chemical composition of the first core corresponds to the initial composition of the 
circumstellar disk; the total of gas and ice abundances at r > 3 AU (T < 500 K) in the 
first core is mostly the same as those in the cold molecular cloud core, while the inner part 
(T > 500 K) is dominated by the simple molecules, such as CO, H 2 and N 2 . 

When the protostar is born, the mass of the circumstellar disk (~0. 01-0.1 Mq) 
overwhelms that of the protostar (~10~ 3 Mq). It causes the gravitational instability, and 
angular momentum is efficiently transported via spiral arms. The mass of the protostar 
and circumstellar disk increase with time via the accretion from the circumstellar disk and 
infalling envelope, respectively. Therefore most of the material constituting the first core 
would accrete to the protostar, and be replaced by the newly accreted envelope material 
during the disk evolution, although a small fraction of the first core material may have been 
transported to outer region of the disk carrying angular momentum. Since the envelope 
material in the protostellar stage is more abundant in large organic molecules than in the 
first core stage (Section 4.3), the circumstellar disk would be enriched with the large organic 
molecule. The large organic molecules also might be formed in the circumstellar disk via 
gas-phase and/or grain-surface reactions if the fluid parcels stay in the disk long enough 
(10 4 -10 5 yrs) (e.g., Charnley et al. 1992; Garrod & Herbst 2006; Aikawa et al. in prep). On 
the other hand, increasing UV radiation and decreasing attenuation in the envelope cause 
the photodissociation of molecu les, such as H 2 Q and CO, at the surfaces of the disk and the 



wall of the outflow cavity (e.g., IVisser et al. 1 120 1 lh . 
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5. CONCLUSION 

We have investigated molecular evolution in the first hydrostatic core stage. We 
performed the three-dimensional radiation hydrodynamic simulation from the molecular 
cloud core to the first core. We traced the fluid parcels, which follow the local flow of 
the gas, and solved molecular evolution along the trajectories, coupling gas-phase and 
grain-surface chemistry. We derived the spatial distributions of molecular abundances and 
column densities in the molecular cloud core harboring the first core. Our conclusions are 
as follows. 

1. We found that the total molecular abundances of gas and ice of many species in 
the cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas 
abundances in the warm envelope and outer layers of the first core (T < 500 K) are mainly 
determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, 
such as H 2 CO, start to be destroyed mainly via collisional dissociation or reactions with H 
atom, and converted to the simple molecules, such as CO, H 2 and N 2 . On the other hand, 
some molecules are effectively formed; carbon-chains, such as C2H2 and cyanopolyynes, are 
formed at the temperature of >700 K. 

2. Gaseous large organic molecules are associated with the first core (r < 10 AU). 
Although the abundances of large organic molecules in the first core stage are comparable 
or smaller than in the protostellar hot corino, their column densities in our model are 
comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. 
We propose that large organic molecules can be good tracers of the first cores. 

3. The shock wave occurs when the envelope material accretes onto the first core. We 
found that the shock heating has little effect on chemistry, if the maximum temperature 
in the shock heating layer is below 1000 K. Even if the gas temperature reaches 2000 K, 
the heating has a limited effect on chemistry. While it significantly decreases the H 2 CO 
abundance and increases the C 2 H 2 abundance in the outer parts of the first core (r < 3 
AU), it has little effect on other species. 

4. Recent hydrodynamic simulations have shown that outer region of the first core 
could directly evolve to the circumstellar disk. In that case, the chemical composition of 
the first core corresponds to the initial composition of the circumstellar disk; the total 
abundances of gas and ice at r > 3 AU (T < 500 K) are mostly determined in the cold 
molecular cloud cores, while inner parts (T > 500 K) are dominated by simple molecules, 
such as CO, H 2 and N 2 . 
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Fig. 1. — Temporal variations of gas number density (solid line) and temperature (dashed 
line) at the center of the core. 
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Fig. 2. — Distributions of gas number density (a, b, d, e) and temperature (c, f ) at £fc = 1500 
yr (a-c) and 2800 yr (d-f). The panels (a) and (d) represent the cross sections at z = 0, 
while the other panels represent the cross sections at y — 0. The arrows denote the velocity 
fields. The white dashed lines in the panels (c) and (f ) depict isothermal lines of 25 K and 
100 K. 
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Fig. 3. — Distributions of gas number density, radial velocity (v r ) and temperature perpen- 
dicular to (a-c) and along the rotational axis (d-f) at tpc — 1500 yr (solid lines) and 2800 
yr (dashed lines). Distributions perpendicular to the rotational axis are in the z = plane 
and averaged in the azimuth direction. 
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Fig. 4. — Spatial distributions of fluid parcels in the middle stage (left) and the late stage 
(right). The parcels are color-coded according to their initial radius (above the equatorial 
plane) and initial specific angular momentum (below the equatorial plane). Gray scale shows 
the distributions of gas number density averaged in the azimuth direction. A parcel is plotted 
at (R, z) = (y/Xp + Up, z p ), when the position of the parcel is (x p , y p , z p ). R < regions are 
simply copied from R > regions. 
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Fig. 5. — Trajectory of the fluid parcel, which is initially at 1300 AU from the center and 
70° degrees from the rotational axis, plotted over the distribution of the gas number density 
when the parcel enters the first core (gray scales and contour lines). The values on the 
contour lines are number density in logarithm scale. The panels (a) and (b) show the whole 
trajectory, while the panels (c) and (d) are the closeup views. The trajectory is projected 
to the z = plane in the panels (a) and (c), and to the y = plane in the panels (b) and 
(d). The trajectory is color-coded according to the temperature in the parcel; the cold phase 
(10-20 K, blue), the warm-up phase (20-200 K, green) and the hot phase (>200 K, red). In 
the panels (c) and (d), the white dashed lines depict the first core surface, and the diamonds 
depict where the parcel passes the accretion shock at the first core surface. 
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Fig. 6. — Temporal variations of the position of the parcel in r, R, z (a), temperature (b, d), 
number density of H nuclei (c), and molecular abundances (e-h) in the fluid parcel. The 
horizontal axis is the logarithm of t^ na x — t. In order to zoom-up the warm-up phase, the 
panels (c-h) are divided into three phases; the cold phase (T < 20 K), the warm-up phase 
(20 K < T < 200 K) and the hot phase (T > 200 K). Note that the scale of the horizontal 
axis is different among the phases. The vertical dotted lines represent the time when the 
parcel passes the accretion shock at the first core surface. In the panels (c-h), the solid 
lines represent gas-phase species, while the dashed lines represent ice-mantle species. HC n N 
indicates total molecular abundances of HC5N, HC7N and HC9N. 
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Fig. 7. — Spatial distributions of fluid parcels, which are color-coded according to abundance 
of gaseous CH 3 OH, in the middle stage (a, b) and the late stage (c, d). The panels (a) and 
(c) represent the cross sections at z = 0, while (b) and (d) represent the R-z plane. In the 
panels (a) and (c), the parcels with \z p \ < 0.5 AU are plotted. The gray scales and the 
contour lines represent the gas number density and temperature distributions, respectively 
in logarithmic scale. 
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Fig. 8. — Abundance of gaseous CH 3 OH as a function of temperature in the middle stage 
(gray square) and the late stage (black cross). 




Fig. 9. — Radial distributions of physical parameters (a, b) and molecular abundances (c-h) 
in the z = plane in the middle stage (a, c, e, g) and the late stage (b, d, f, h). In the panels 
(c-h), the solid lines represent gas-phase species, while the dashed lines represent ice- mantle 
species. HC n N indicates total molecular abundances of HC5N, HC7N and HC9N. 




Fig. 10. — Distributions of physical parameters (a, b) and molecular abundances (c-h) along 
the rotational axis. Other details are the same as in Figure [9j 
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Fig. 11. — Column densities along the rotational axis integrated from the outer edge of the 
envelope to z, N z ^{z) = J^°riidz', in the middle stage (upper panel) and the late stage 
(lower panel). The vertical solid lines represent the first core surface. The vertical dashed 
and dotted lines represent the heights from the midplane, where the optical depths of dust 
continuum reach unity with the wave length of 100 //m and 1 mm, respectively. 
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Fig. 12. — Gaseous molecular column densities in front of the Ti mm = 1 plane in the face-on 
object. 
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Fig. 13. — Temperature (upper panel) and total (gas and ice) molecular abundances of 
selected species (lower panel) along the trajectory discussed in Section 3.3 after the parcel 
enter the first core at r ~ 7 AU. The solid and dashed lines in the lower panel represent the 
abundances with and without the shock heating, respectively. 
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Table 1: Generalized reaction scheme of gas-grain interactions. X, Y, and Z are neutral 
gaseous species whose charges are denoted by their indices, and gr are grain particles. X ice 
represents the species on grain surfaces. 
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Table 2: Initial abundances relative to hydrogen nuclei. 



Species Abundance 1 

H 2 5.00(-l) 

H 5.00(-5) 

He 9.75(-2) 

C + 7.86(-5) 

N 2.47(-5) 

1.80(-4) 

Si+ 2.74(-9) 

S+ 9.14(-8) 

Fe+ 2.74(-9) 

Na+ 2.25(-9) 

Mg+ 1.09(-8) 

C1+ 1.00(-9) 

P+ 2.16(-10) 



1 a(— b) means a x 10 
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Table 3: Temperature, number density of hydrogen nuclei, position in r, and duration of 
the cold, warm-up, and hot phases in the fluid parcel, tff is the free-fall time of the initial 
central density. 





Temperature 


?i H 


Position in r 


Duration 




(K) 


(cm -3 ) 


(AU) 


(y r ) 


Cold phase 


10-20 


1 x 10 6 -2 x 10 9 


1500-40 


7 x 10 4 (~3t//) 


Warm-up phase 


20-200 


2 x 10 9 -8 x 10 12 


40-5 


120 


Hot phase 


200-2000 


8 x 10 12 -8 x 10 15 


5-<0.1 


1900 



Table 4: Comparison between the protostellar and prestellar hot corino models. 





Protostellar hot corino 


Prestellar hot corino 


Radius (T = 100 K) 


120 AU a 


<10 AU 


Gas density (T = 100 K) 


10 7 cm" 33 


10 10 cm -3 


H 2 CO 


2.5(-5) 


2.6(-5) 


CH 3 OH 


1.0(-5) 


6.4(-6) 


HCOOCH3 


3.3(-9) 


4.8(-ll) 


CH 3 OCH 3 


6.0(-ll) 


2.2(-13) 


HCOOH 


3.0(-9) 


1.3(-10) 



Note. — As the protostellar HC model, we adopt the values at the protostellar age of 9.3 x 10 4 yr. As the 
prestellar HC model, we use the values in the middle stage. Model abundances are represented with respect 
to H2 molecule at 15 AU in the protostellar HC and at 5 AU in the prestellar HC. 
"Masunaga & Inutsuka (2000). 
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Appendix A : Three Body Association and Collisional Dissociation 

Reactions 

Table Al lists the three body association reactions, while Table A2 lists the collisional 
dissociation reactions in our network. 
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Table Al: Three body association reactions in the reaction network, which are taken from Willacy et al 
(1998). The rate is calculated as k = a(T/300) /3 exp(- 7 /T). 









Reaction 






a 


13 


7 


CH 3 - 


- H - 


h M 




CH 4 


h M 


9.62(-31) 


-1.80 


3.23(3 


C 


O 


h M 


— > 


CO 


M 


2.14(-29) 


-3.08 


2.11(3 


CO 


O 


h M 


— > 


CO2 - 


M 


9.56(-34) 


0.00 


1.06(3 


H 


- CO - 


h M 


— > 


HCO - 


M 


5.30(-34) 


0.00 


3.70(2 


N 


N 


h M 




N 2 


M 


9.53(-34) 


-0.50 


0.00(0 


N 


H 


h M 




NH 


M 


3.09(-32) 


-0.50 


0.00(0 


NH 


H 


h M 




NH 2 


M 


3.18(-33) 


-0.50 


0.00(0 


NH 2 - 


H 


h M 




NH 3 


M 


1.32(-33) 


0.00 


4.37(3 


N 


O 


h M 




NO 


M 


3.86(-34) 


0.00 


2.57(3 


O 


O 


h M 




O2 


M 


5.25(-35) 


0.00 


9.02(2 


O 


H 


h M 




OH 


M 


6.32(-34) 


0.00 


1.96(3 


OH 


H 


h M 




H 2 - 


h M 


6.86(-31) 


-2.00 


0.00(0 


s 


S 


h M 


— > 


s 2 


h M 


2.76(-33) 


-0.50 


0.00(0 


so 


O 


h M 


— > 


so 2 - 


h M 


1.86(-31) 


-0.50 


0.00(0 


H 


H 


h H 2 




H 2 


- H 2 


8.72(-33) 


-0.60 


0.00(0 


H 


H 


h H 




H 2 


- H 


1.83(-31) 


-1.00 


0.00(0 



-49 - 



Table A2: Collisional dissociation reactions in the reaction network. The rate is calculated as k = 
a(T/300)^ exp(— 7/T). At the end of the row, W and H indicate that the rate is taken from Willacy et 
al. (1998) and Harada et al. (2010), respectively. WE also refers Willacy et al., but 7 is replaced by the 
difference between the formation enthalpy of the products and the reactants at 300 K. E indicates that the 
rate is estimated. The rates of the collisional dissociation reactions, which are not listed here, are calculated 
as k = 10~ 10 exp(-50000/r). 



Reaction a /3 7 Reference 



c 2 


(- M 


— > 


C 


+ 


C 


+ M 


6.17(-10) 


0.00 


7.17(4) 


WE 


c 3 


1- M 


-> 


c 2 


+ 


C 


+ M 


6.17(-10) 


0.00 


8.83(4) 


WE 


c 4 


1- M 


-> 


c 3 


+ 


C 


+ M 


l.OO(-lO) 


0.00 


6.81(4) 


E 


c 5 


1- M 


— V 


c 4 


+ 


C 


+ M 


l.OO(-lO) 


0.00 


8.52(4) 


E 


CH 


1- M 


— > 


c 


+ 


H 


+ M 


3.16(-10) 


0.00 


4.10(4) 


WE 


CH 2 


1- M 


— > 


CH 


+ 


H 


+ M 


6.64(-9) 


0.00 


5.12(4) 


WE 


CH 2 


I- M 


— ► 


c 


+ 


H 2 


+ M 


2.16(-10) 


0.00 


3.97(4) 


WE 


CH 3 


1- M 


->• 


CH 


+ 


H 2 


+ M 


1.15(-9) 


0.00 


5.40(4) 


WE 


CH 3 


I- M 


— > 


CH 2 


+ 


H 


+ M 


1.70(-8) 


0.00 


5.52(4) 


WE 


CH 4 


I- M 


— ► 


CH 3 


+ 


H 


+ M 


1.20(-6) 


0.00 


5.28(4) 


WE 


C 2 H 


1- M 


— > 


2 


+ 


H 


+ M 


3.74(-10) 


0.00 


6.96(4) 


WE 


C 2 H 2 


1- M 


— > 


C 2 


+ 


H_> 


+ M 


7.59(-17) 


0.00 


7.35(4) 


WE 


C 2 H 2 


I- M 


— > 


C 2 H 


+ 


H 


+ M 


6.97(-8) 


0.00 


5.63(4) 


WE 


C 2 H 3 - 


1- M 


— ► 


C 2 H 2 


+ 


H 


+ M 


5.07(-20) 


-7.50 


2.29(4) 


W 


C 2 H 4 


1- M 


->• 


C 2 H 3 


+ 


H 


+ M 


l.OO(-lO) 


0.00 


5.59(4) 


E 


C 2 H5 


I- M 


— ► 


C 2 H 4 


+ 


H 


+ M 


l.OO(-lO) 


0.00 


1.67(4) 


E 


C 3 H 4 - 


1- M 


— > 


C 3 H 3 


+ 


H 


+ M 


1.66(-7) 


0.00 


4.47(4) 


WE 


C 4 H 2 


1- M 


— »• 


C 4 H 


+ 


H 


+ M 


5.91 (-7) 


0.00 


4.03(4) 


W 


CO 


I- M 


— >• 


C 


+ 





+ M 


1.93(-6) 


0.00 


1.29(5) 


WE 


co 2 


I- M 


— > 


CO 


+ 





+ M 


8.01(-11) 


0.00 


6.40(4) 


WE 


HCO 


1- M 


— ► 


CO 


+ 


H 


+ M 


2.39(-10) 


0.00 


8.13(3) 


W 


H 2 CO - 


I- M 


— > 


HCO 


+ 


H 


+ M 


2.10(-8) 


0.00 


4.54(4) 


WE 


H 2 CO - 


I- M 


— ► 


H 2 


+ 


CO 


+ M 


3.95(-9) 


0.00 


1.75(4) 


W 


c 2 o 


1- M 


— > 


C 


+ 


CO 


+ M 


l.OO(-lO) 


0.00 


3.84(4) 


E 


CS 


I- M 


— >• 


C 


+ 


s 


+ M 


l.OO(-lO) 


0.00 


8.58(4) 


E 


HCS 


I- M 


— ► 


CS 


+ 


H 


+ M 


l.OO(-lO) 


0.00 


2.39(4) 


E 


H 2 CS 


1- M 


->• 


HCS 


+ 


H 


+ M 


l.OO(-lO) 


0.00 


4.85(4) 


E 


CN 


1- M 


— »• 


C 


+ 


N 


+ M 


3.32(-10) 


0.00 


9.07(4) 


WE 


HCN 


I- M 


— ► 


CN 


+ 


H 


+ M 


2.08(-8) 


0.00 


6.23(4) 


WE 


HNC 


1- M 


— > 


CN 


+ 


H 


+ M 


l.OO(-lO) 


0.00 


5.52(4) 


E 


HNCO - 


1- M 


— > 


NH 


+ 


CO 


+ M 


l.OO(-lO) 


0.00 


4.42(4) 


E 


C 2 N 


I- M 


— >• 


C 


+ 


CN 


+ M 


l.OO(-lO) 


0.00 


7.16(4) 


E 


S 2 


I- M 


— ► 


s 


+ 


S 


+ M 


7.95(-ll) 


0.00 


5.12(4) 


WE 


HS 


1- M 




11 


+ 


s 


+ M 


l.OO(-lO) 


0.00 


4.28(4) 


E 


H 2 S 


I- M 


— > 


11s 


+ 


H 


+ M 


7.70(-10) 


0.00 


4.55(4) 


WE 


SO 


I- M 


— ► 


S 


+ 





+ M 


6.61(-10) 


0.00 


6.27(4) 


WE 


so 2 


1- M 


— > 


so 


+ 





+ M 


6.60(-9) 


0.00 


6.63(4) 


WE 


H 2 S 2 


I- M 


— > 


H 2 S 


+ 


s 


+ M 


l.OO(-lO) 


0.00 


2.90(4) 


E 


N 2 


I- M 


— ► 


N 


+ 


N 


+ M 


9.21(-5) 


-2.50 


1.14(5) 


WE 


NH 


1- M 




N 


+ 


H 


+ M 


4.40(-10) 


0.00 


3.80(4) 


W 


NH 2 


1- M 




NH 


+ 


H 


+ M 


1.00(-36) 


0.50 


4.91(4) 


WE 


NH 3 


I- M 


— ► 


NH 2 


+ 


H 


+ M 


3.46(-8) 


0.00 


5.42(4) 


WE 


NH 3 


1- M 


— > 


NH 


+ 


H 2 


+ M 


1.05(-9) 


0.00 


5.08(4) 


WE 


NS 


1- M 


— > 


N 


+ 


S 


+ M 


l.OO(-lO) 


0.00 


5.85(4) 


E 


NO 


I- M 


— ► 


N 


+ 





+ M 


4.10(-9) 


0.00 


7.60(4) 


WE 


HNO 


1- M 




11 


+ 


NO 


+ M 


l.OO(-lO) 


0.00 


2.36(4) 


E 


N 2 


1- M 


->• 


N 


+ 


NO 


+ M 


l.OO(-lO) 


0.00 


5.79(4) 


E 


SiH 


I- M 


— > 


Si 


+ 


H 


+ M 


l.OO(-lO) 


0.00 


3.51(4) 


E 


SiH 2 


I- M 


— ► 


SiH 


+ 


H 


+ M 


7.60(-8) 


-1.80 


3.87(4) 


WE 


SiH 2 


1- M 


— > 


Si 


+ 


H_> 


+ M 


7.60(-8) 


-1.76 


2.13(4) 


WE 


SiH 3 


I- M 


— ► 


SiH 


+ 


H 2 


+ M 


4.98(-10) 


0.00 


2.12(4) 


WE 


SiH 4 


I- M 


— ► 


SiH 3 


+ 


H 


+ M 


5.00(-ll) 


0.00 


4.62(4) 


WE 


SiH 4 


1- M 


— »• 


SiH 2 


+ 


H 2 


+ M 


5.00(-ll) 


0.00 


2.90(4) 


WE 
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Table A2: continue 







Reaction 






a 


P 


7 


Reference 


SiC 


h M 


-> Si 


C 


h M 


1.00(-10) 


0.00 


5.38(4) 


E 


SiS 


- M 


-> Si 


s 


- M 


1.00(-10) 


0.00 


7.47(4) 


E 


SiN 


h M 


-> Si - 


N 


h M 


1.00(-10) 


0.00 


6.62(4) 


E 


SiO 


- M 


-> Si 


O 


- M 


1.00(-10) 


0.00 


9.62(4) 


E 


Si0 2 


- M 


-> SIO 


O 


- M 


1.00(-10) 


0.00 


5.46(4) 


E 


SiC 2 


h M 


-> SiC 


C 


h M 


1.00(-10) 


0.00 


9.95(4) 


E 


OH 


- M 


-> 


H 


- M 


4.00(-9) 


0.00 


5.15(4) 


WE 


H 2 


- M 


-> OH 


H 


- M 


5.80(-9) 


0.00 


6.00(4) 


WE 


OCS 


h M 


-> CO 


S 


h M 


4.82(-10) 


0.00 


3.71(4) 


WE 


HC1 


- M 


-> H 


CL 


- M 


1.00(-10) 


0.00 


5.19(4) 


E 


CC1 


- M 


-> c - 


CI 


- M 


1.00(-10) 


0.00 


4.04(4) 


E 


CIO 


h M 


-> CI - 


O 


h M 


1.00(-10) 


0.00 


3.24(4) 


E 


CP 


- M 


-> c - 


P 


- M 


1.00(-10) 


0.00 


7.02(4) 


E 


HCP 


- M 


-> H 


CP 


- M 


1.00(-10) 


0.00 


6.23(4) 


E 


CH 2 PH 


- M 


-> CH 2 


PH 


- M 


1.00(-10) 


0.00 


1.11(5) 


E 


MgH 


- M 


-> Mg 


H 


- M 


1.00(-10) 


0.00 


2.36(4) 


E 


NaH 


- M 


-»• Na 


H 


- M 


1.00(-10) 


0.00 


2.42(4) 


E 


PH 


h M 


-> P 


H 


h M 


1.00(-10) 


0.00 


3.38(4) 


E 


PH 2 


h M 


-> PH 


H 


h M 


1.00(-10) 


0.00 


4.13(4) 


E 


PN 


h M 


-> P 


N 


h M 


1.00(-10) 


0.00 


8.23(4) 


E 


PO 


h M 


-> P 


O 


h M 


1.00(-10) 


0.00 


7.09(4) 


E 


HPO 


h M 


->■ H 


PO 


h M 


1.00(-10) 


0.00 


3.02(4) 


E 


o 2 


h M 


-► O 


O 


h M 


5.16(-10) 


0.00 


6.00(4) 


WE 


2 H 


h M 


-> o 2 


H 


h M 


1.00(-10) 


0.00 


2.47(4) 


E 


H 2 2 


h M 


-+ H 2 - 


O 


h M 


1.00(-10) 


0.00 


1.73(4) 


E 


O3 


h M 


-> o 2 


O 


h M 


1.00(-10) 


0.00 


1.28(4) 


E 


CH 3 N 


h M 


-> CH 3 - 


N 


h M 


1.00(-10) 


0.00 


6.34(4) 


E 


CH 3 NH 2 - 


h M 


-> CH 3 - 


- NH 2 - 


h M 


1.00(-10) 


0.00 


4.26(4) 


E 


CH 2 CN 


h M 


-> C 2 N - 


H 2 


h M 


1.00(-10) 


0.00 


6.00(4) 


E 


CH3CN 


h M 


->■ CH 3 - 


CN 


h M 


1.00(-10) 


0.00 


6.21(4) 


E 


HCOOH 


h M 


-> HCO - 


OH 


h M 


1.00(-10) 


0.00 


5.55(4) 


E 


CH 2 CO 


h M 


->■ CH 2 - 


CO 


h M 


5.97(-9) 


0.00 


3.86(4) 


WE 


CH3OH 


h M 


->■ CH 3 - 


OH 


h M 


3.32(-7) 


0.00 


4.64(4) 


WE 


CH3CHO - 


h M 


-> CH 3 - 


- HCO - 


h M 


1.00(-10) 


0.00 


4.20(4) 


E 


H 2 


h e- 


->■ H 


H 


h e - 


2.00(-9) 


0.50 


1.20(5) 


H 


H 2 


h H 


->■ H 


H 


h H 


1.00(-10) 


0.00 


5.20(4) 


H 


H 2 


h He 


-)■ H 


H 


h He 


l.OO(-ll) 


0.00 


5.20(4) 


11 


H 2 


h H 2 


-)■ H 


H 


h H 2 


1.30(-11) 


0.00 


5.20(4) 


H 



